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ABSTRACT 


This implementation of the physical optics approximation method for predicting 
the radar cross section (RCS) of complex objects utilizes the scientific computational 
features of MATLAB and its Graphical User Interface (GUI) functions to provide an 
error-free encoding of input parameters and efficient calculation. The POFACETS GUI 
provides a convenient tool for both students and electromagnetic professionals to design 
complex models by representing its component parts as triangular facets and offers 
options for defining the surface roughness of the scattering object. The GUI calculates 
the monostatic or bistatic radar cross section of the modeled object based on certain radar 
frequencies and parameters of interest, creates a variety of plots for the model geometry 
and its RCS measurement, and provides additional functionalities to enhance 
visualization of the object and verify correctness of data. Various help windows and 
messages are available to the user to serve as guide in the efficient use of the GUI. The 
simplicity of the GUI design does not require the user to have a detailed understanding of 
the various steps involved in the RCS calculation. 
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I. INTRODUCTION 


Designing operational platforms and predicting their radar cross section (RCS) for 
a range of frequencies provides a challenging problem in modem electronic warfare. One 
approach to high frequency prediction calculations is to estimate a complex model with 
an array of simple shapes, such as triangles or flat plates. The RCS is obtained by 
computing the scattered field of the collection of these simple shapes to obtain the total 
RCS of the target. 

Among the most common RCS prediction methods for any arbitrary three- 
dimensional target are the Method of Moments (MOM), the Finite Difference Method 
(FDM), Geometrical Optics (GO) and Physical Optics (PO). The MOM is a popular 
technique that solves the electric or magnetic field integral equation by reducing it to a 
matrix problem with the matrix size directly related to the electrical size of the scattering 
body in terms of the incident wave. This method gives very accurate results but is 
computationally intense. In the FDM, finite differences are used to approximate the 
differential operators in Maxwell’s equations in either the time or frequency domain. The 
target must be discretized as in the MOM. Since this method necessarily calculates the 
fields in the region around the target, computing the RCS of a target with a characteristic 
dimension of several orders of magnitude of the wavelength would entail considerable 
time. GO and the Geometric Theory of Diffraction (GTD) involve ray tracing methods to 
analyze electrically large targets of arbitrary shape. Diffracted rays are allowed to 
originate from the scattering of the incident wave at edges, comers, and vertices. Among 
the problems associated with GO is that the source and the observation point must satisfy 
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Snell’s Law or the reflected field is zero. It yields discontinuous fields at shadow 
boundaries and zero fields in the shadow regions. GTD solves this problem by 
accounting for the edge-related effects. PO overcomes the problems associated with the 
infinities of ray tracing for flat and singly curved surfaces by approximating the induced 
surface fields. The induced currents are integrated over the illuminated portions of the 
target to obtain the scattered far field while setting the current to zero over the shadowed 
portions. It is most accurate in the specular direction and gives good estimates at points 
far from discontinuities and shadow boundaries. However, because PO abruptly sets the 
current to zero at a shadow boundary, the computed field values at wide angles and in the 
shadow regions are inaccurate. Also, surface waves, multiple reflections and edge 
diffractions are not included [Refs. 1,2 and 3]. 

There is a need for an easy-to-use tool that conveniently allows the modeling of 
an arbitrary three-dimensional object with component triangular facets, visualization of 
the geometry, and calculation of the radar cross section. The physical optics 
approximation is an efficient method that gives accurate results for most large targets, yet 
is not too computationally demanding. To make this tool appealing to the user and 
prevent, if not eliminate unnecessary errors in the setting up of various parameters, the 
design of a graphical user interface is desirable. 

This thesis focuses on the design, development, and implementation of a 
Graphical User Interface (GUI) in MATLAB for a PO-based RCS prediction code. The 
program allows creation of a model comprised of triangular facets with options for 
specifying the surface characteristics. It then calculates the monostatic or bistatic radar 
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cross sections given certain computational parameters without requiring the user to have 
a detailed understanding of the various steps involved in the calculation. 

The end product of this thesis is the POFACETS GUI that implements a physical 
optics approximation method that may be used by both students and electromagnetic 
professionals as a convenient tool to design complex models and predict their RCS 
performance based on the radar frequency and other parameters of interest. The GUI 
provides additional functionalities to enhance visualization of the model geometry and 
generates a variety of plots to present the RCS data. The GUI is based on codes 
developed by Professor David Jenn and utilizes the scientific computational features of 
MATLAB. 

Chapter II gives a brief description of the theory behind the formulation of the 
code by providing an overview of radar fundamentals, electromagnetic scattering and a 
comparative discussion of the common methods of high frequency RCS prediction. 

In Chapter in the flow or framework of the code is presented with the supporting 
mathematics involved in using the physical optics method to calculate the RCS. 

Chapter IV describes the POFACETS Graphical User Interface. It is intended to 
provide a descriptive guide in making the user acquainted with the operation of the GUI. 
A step-by-step illustration is provided for working on the examples that are discussed. 

Finally, Chapter V summarizes the research and suggests further work that could 
be performed to enhance the code. 
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II. THEORETICAL FORMULATION 


A. RADAR BACKGROUND 

The acronym radar was coined during World War II from the phrase radio 
detection and ranging. A radar is a device that transmits an electromagnetic wave and 
detects objects by virtue of the energy scattered from them in the direction of the 
receiver. Generally, the radar transmitter and receiver can be at different locations when 
viewed from the target and this is called a bistatic radar. In most applications the 
transmitter and receiver are co-located and often share the same antenna, a configuration 
referred to as a monostatic radar. The classical form of the radar equation is 


P r = (•£%)(— 
4nR 47tR 


^ A rrr 


where P r is the power received (in watts), P t the output power of the transmitter, G t the 
gain of the transmitter, G r the gain of the receiver, crthe radar cross section (in square 


meters), A the wavelength (in meters), L the reduction in the received signal due to losses 
(0 < L < 1), and R the range of the target from the radar. The first term in parenthesis is 
the power density at the target (in watts/m 2 ) due to radiation from the transmitter. The 
second term is the incident power reflected or scattered back to the receiver. The third 
term is the amount of the reflected power which is captured by the receiving antenna due 
to its effective area. For the case of a monostatic radar with the same antenna used for 


transmitting and receiving (G ( is the same as G r ), equation 2.1 simplifies to 

P&oXh 
r ~ (4 7lfR* 


( 2 . 2 ) 
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B. ELECTROMAGNETIC SCATTERING 


1. Radar Cross Section 

Radar cross section (RCS) is a measure of the power scattered in a given direction 
when a target is illuminated by an incident wave, normalized to the power density of the 
incident field. The purpose of the normalization is to remove the effect of the range and 
arrive at an RCS computation that is independent of the distance between the target and 

illumination source. In terms of the incident and scattered electric field intensities, E t 

and E s , the RCS is defined as 


cr= ii m 4 nR 2 

R —)■ 



(2.3) 


I *1 

The unit of RCS most commonly used is decibels relative to a square meter (dBsm): 

cr, dBsm = 101og(<j,m 2 ) (2.4) 

Radar cross section is a function of many factors which include the target configuration 
and its material composition, frequency or wavelength, transmitter and receiver 
polarizations, and the target aspect (angular orientation of the target) relative to the radar. 
Hence, in general, a could be specified as a pq (6,0) where q and p pertain to the incident 
and received polarizations, respectively, and where (6,0) are the polar spherical angles of 
view as shown in Figure 1. 
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Figure 1. Coordinate System for Incident Polarization 


2. Scattering Regions 

The scattering characteristics of a target are dependent on the frequency of the 
incident wave. Three regimes distinctly differentiate the RCS scattering of a target 
depending on the ratio of the wavelength X to the body size L. The first is the low 
frequency region or Rayleigh region (A » L) where the phase variation of the incident 
plane wave across the target is small and the shape of the body is not of importance. The 
second is the resonance region or Mie region (A — L) where the phase variation of the 
current across the body is significant and all parts contribute to the scattering pattern. The 
third is the high frequency region or Optical region (A « L) where the phase of the 
current across the body varies with many cycles and the scattering process can often be 
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considered as the summation of the returns from non-interacting isolated scattering 
centers [Refs. 1,2 and 4], 

3. Target Scattering Matrix 

The radar cross section is a scalar number that is a function of the polarization of 
the incident and received waves. The polarization of the incident wave from the radar 
will not always be in the same orientation as the polarization of the target coordinate 
system and thus, it is necessary to decompose the polarization of the incident wave into 
components in the target coordinate system. Since the wave vector components at large 
distances are tangent to a sphere, two orthogonal components in terms of the variables 6 
and <f> are sufficient to represent the incident field in a spherical system centered at the 
target. Thus, in general, the incident field can be written as 

E, = E w d + ( 2 . 5 ) 

A A 

where 6 and 0 are unit vectors in the target coordinate system. The polarization of the 
scattered field is not necessarily the same as the incident field since most complex targets 
generate a cross-polarized scattering component due to multiple diffractions and 
reflections. The following scattering matrix can be used to specify the relationship 
between the incident and scattered fields: 

E s e See^etp E ig 

E S S E (2*6) 

s<j> J _ 00 ^ 00 ^i 0 

The S pq are the scattering parameters, where the first index specifies the 

polarization of the receive antenna and the second, the polarization of the incident wave. 
The elements of the scattering matrix are complex quantities and in terms of the RCS 
[Ref. 1], 
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(2.7) 


S 


pq 


JEH j<r* 

y[4ftR 2 




= tan 


' ¥W ' 

^Re(£ jp /£,,)“ 


B. COMMON RCS PREDICTION METHODS 

Among the most common numerical RCS prediction methods for any arbitrary 
three-dimensional target are the Method of Moments (MOM), the Finite Difference 
Method (FDM), Geometrical Optics (GO) and Physical Optics (PO). The first two 
provide rigorous solutions but are computationally intense requiring considerable 
computer time to calculate results, while the latter two are approximate methods which 
give accurate and fast results, but have limitations. 

1. Method of Moments 

The MOM involves breaking the target into a number of surface or line 
subdomains. This popular solution technique solves the electric or magnetic field 
integral equation by reducing it to a matrix problem with the matrix size directly related 
to the electrical size of the scattering body in terms of the incident wavelength. The 
coefficients of a series approximation to the current can then be determined using matrix 
algebra. In general, the subdomains must be small compared to the wavelength. This 
method gives very accurate solutions and is used in programs such as the Numerical 
Electromagnetics Code (NEC) [Ref. 3]. However, the calculation time is dictated by that 
required to solve the simultaneous equations and the method may not be practical for 
large targets at high frequencies [Refs. 1 and 3]. 
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2 . 


Finite Difference Methods 


In this method, finite differences are used to approximate the differential 
operators in Maxwell’s equations, usually in the time domain. The target must be 
discretized as in the MOM and the elements have to be small compared to the 
wavelength. Finite difference techniques do not require the solution of large dense 
matrices that the MOM does because the solution is stepped in time throughout the 
scattering body. Since this method calculates the fields in a computational grid around 
the target, calculating the RCS of a target with a characteristic dimension of several 
orders of magnitude of the wavelength would entail considerable time to compute [Refs. 
1 and 3]. 

3. Geometrical Optics 

GO involves ray tracing methods to analyze electrically large targets of arbitrary 
shape. The rules for ray tracing in a simple medium (linear, homogeneous, and isotropic) 
are similar to reflection and refraction in optics. In addition, diffracted rays are allowed 
that originate from the scattering of the incident wave at edges, comers, and vertices. 
The formulas are derived on the basis of infinite frequency (A —> 0) which implies an 
electrically large target. Among the problems associated with GO is that the source and 
the observation points must satisfy Snell’s Law or the reflected field is zero. It yields 
discontinuous fields at shadow boundaries and zero fields in the shadow regions. The 
Geometric Theory of Diffraction was developed to take into account the edge-related 
effects by decomposing the total field at the observation point into its GO and diffracted 
components, thereby resulting in approximations for edge diffracted fields [Refs. 1 and 
4]. 
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4. Physical Optics 

Physical Optics overcomes the problems associated with the infinities of ray 
tracing for flat and singly curved surfaces by approximating the induced surface current 
on the illuminated portions of the body to be proportional to the incident magnetic field 
intensity and integrating them to obtain the scattered far field. On the shadowed portion 
of the target, the current is set to zero. Thus, the current is approximated as 

[in x H i , for the illuminated portion 
J s — \ (A®) 

[0, for the shadowed portion 

where H i is the incident magnetic field intensity at the surface and n is the surface 
normal unit vector. 

The approximation in equation (2.8) is used in the radiation integrals to compute 
the scattered far field from the target. PO is a high-frequency approximation method that 
gives best results for electrically large bodies (L > 1C )A ). It is most accurate in the 
specular direction and gives good estimates at points far from discontinuities and shadow 
boundaries. However, because PO abruptly sets the current to zero at a shadow 
boundary, the computed field values at wide angles and in the shadow regions are 
inaccurate. Also, surface waves, multiple reflections and edge diffractions are not 
included [Refs. 1, 2 and 4]. 
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III. PHYSICAL OPTICS APPLIED TO TRIANGULAR PATCHES 


A. SCATTERING FROM FACETS WITH ARBITRARY ORIENTATION 

The geometry of a single facet with an arbitrary orientation is shown in Figure 2. 
The facet vertices are known in terms of their cartesian coordinates ( x n ,y n ,z „) for 

n = 1...3. It is assumed that the vertices are labeled in a right-hand sense 
(counterclockwise) so that the thumb of the right hand gives the outward normal 
direction. For a closed body that surrounds the origin, the outward normal will in general 
point away from the origin. The position vectors to the nodes are given by 
r n =x n x + y n y + z n Z. 




Figure 2. Arbitrarily Oriented Facet, (a) Arbitrarily Oriented Triangular Facet, 

(b) Facet Sub-Areas for Natural Coordinates 

The edge vectors are determined from the node coordinates by 

\=r3-r 2 , C 3 - 1 ) 
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The tip of edge vector n originates at node n. The outward normal is obtained by taking 
the cross product of any two edge vectors in a right-handed sense. For example, 

A _ 

n-jg = n x x+n y y + n z z (3.2) 

With this ordering of the vertices, the outward side of the facet is denoted as the front 
face and the opposite side as the back face. If a plane wave is incident from an angle 
propagating towards the origin, then its propagation vector is given by 

A t * 

k t =~r = ~(x sin Q l cos fa + y sin 6 X sin fa + z cos 0 t ) = -{xu t + yv ; - + zw t ) (3.3) 

where (m, , v £ -, w t ) are the x, y and z direction cosines, respectively, and r the radial unit 
vector from the origin to the source at (0 iy fa). It is a simple matter to determine 
whether the front face is illuminated or not by performing the test 

-ki»n> 0 (3.4) 

The total area of the facet, A, can be computed from a vector relationship between edges 

A = |pl x ^| (3.5) 

In order to simplify integration over the surface of a triangle, normalized area 
coordinates (also known as simplex or natural coordinates) are used. Consider an 
integration point P located on the plate at (x p ,y p ,Zp) . Its position vector is 

r p =x P x+y P y + z P z (3.6) 

and 

P\=r p -r[, p 2 =r P - r 2 , P 2 =r p - r 3 (3.7) 
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Note that for the monostatic case 


g=h 

The scattered field is written as 


Es 2 nr 


jk 


jkr 


V^e-Ej}\ \e^ds f 


(3.15) 


(3.16) 


POLARIZATION A 
FACTOR, PF '- 




It is not possible to obtain an exact closed form result for the integral. However, 
given that the incident wavefront is plane and known at the facet vertices, the amplitude 
and phase at interior integration points can be found by interpolation. Then the integrand 
can be expanded using a Taylor series, and each term integrated to give a closed form 
result. Usually a small number of terms (on the order of 5) will give a sufficiently 
accurate approximation for the integral. 

Using the notation of [Ref. 5], define 


h =jjc c (T J ,Z)e jD c^ ) dr } dZ 


where 


C c (JJ, €) = C p r) + C q £ + C Q = 
For a unit amplitude plane wave 


Ei 


Similarly, 


E t -1—» C p -C q -0, C Q =1 


D c (tj,£) = D 7]+ D q £ + D 0 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


where 


D p = *[(*! -x 3 )u + (yj -y 3 )v + (zi -z 3 )w] 
D q = ^[(^2 -X 3 )u + (y 2 -y 3 )v + (z2- Z 3 )w] 
D a = k[x 3 u + y 3 v + z 3 w] 

16 


(3.21) 



Now the integral is given by 


\j d p 

C o 

_J D q 

Co 

C 0 

1 

D p {D q -D p )_ 

£ 

D q (D q -Dp)_ 

D q D p\ 


(3.22) 


To avoid the numerical problems that occur near the singularities (denominators near 
zero), the following special cases are used: 

Case 1:1 D, 


h = 


<L t and |D ? | >L t 
2 Ae jD ° ~ 

. _ 2Lt 


0 Dp) 


f -Co 

n\ 

[tt + 1 


+ e J q C 0 G(n,-D a ) 


Case 2: 


D, 


< L t and D q 


< L, 


l e =2Ae^l 


C 0 (jD p Y{jD„Y’ 


Case 3: 


D, 


n =0m=0 (m + n + 2)\ 

>Lf and D q <L t 

jD.JD, ” () D qT( C„ 


l c =2Ae^e' " % 


n=0 


n\ [ n +1 


G(n + 1,—D p ) 


Case 4: \D. 


4 = 


— Lf, D q > L t , and D q Dp < L t 


2Ae»o ~ O' D p -jD q f 


jDq n= 0 ra! 

The function G is defined by 

1 

G(n,y) = j s n e^ s ds 


e iDq C n 


-C 0 G(n,D q ) + -- 

^ n + 1 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


0 

There is recursive relationship that is efficient for finding G for higher values of n 

n>\ 

jy 


(3.28) 
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and 


e jr -1 
G(0,y) = - - t 

jy 


(3.29) 


C. COORDINATE TRANSFORMATIONS 

In general, the local facet coordinate system will not be aligned with the global 
coordinate system. A series of transformations must be performed in order to obtain the 
required relationships between the local and global coordinate variables and unit vectors. 
Figure 3 shows a facet with arbitrary orientation. The global coordinates are (x, y, z) and 
the local coordinates are (x",y",z "). Figure 4 shows that the z" axis is aligned with the 
facet normal. Also, the figure shows a right triangle facet with the x" and y" axes 
aligned with the triangle edges, however, the transformation can be applied to any 
arbitrary facet that lies in the x" - y" plane. 



Figure 3. Facet with Arbitrary Orientation 
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Figure 4. Facet Rotation Angles 


The double-primed variables can be obtained by a series of two rotations of the 
global coordinate variables. The first is a rotation about the z axis by an angle a as 
shown in Figure 5. The resulting coordinate axes are referred to as (x\y',z'). The 
rotation angle is a- atm2(n y ,n x ) (atan2 is the MATLAB arctangent function that gives 

the proper quadrant). The second rotation is by an angle /? = cos _1 (£ • n) about the y' 
axis, illustrated in Figure 6. 


/ 



Figure 5. First Rotation by an Angle a about the z axis. 
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Figure 6. Second Rotation by an Angle (3 about the y' axis. 


The transformation matrices are determined from geometry: 


r /“ 

X 


cos a 

sin a 

O' 

X 

/ 

y 

= 

-sin or 

cos a 

0 

y 

/ 

_Z _ 


0 

0 

1 

_z_ 


T' 


r 

X 


cos/? 

0 

-sin /3 

- /' 
X 

y 

— 

0 

1 

0 

/ 

y 

rt 

Z _ 


sin ft 

0 

cos P 

/ 

_Z _ 




The double overbar denotes a matrix. Therefore, 


|V1 


X 

i- 

M 

* * 
l_ 

_ rr 

y 

_z_ 


The same transformation matrix holds for the unit vectors 


V 

X 


" A “ 

X 

* rr 

y * 

=rr 

y 

X. 


i_ 


(3.30) 


(3.31) 


(3.32) 


(3.33) 
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and the direction cosines 


|V1 


u 

I 

* V 

-1 

=rr 

V 

w 


The direction cosines are 


u = sin # cos (p 
v = sin (9 sin ^ 
w = cos# 


(3.34) 


(3.35) 


and 

u = sm0 cos (p 

v' 1 ' = sin#'sin^' , (3.36) 

w / ” = cos#' r 

These quantities can be either incident angles (subscript “z”) or observation angles 
(subscript “/?”)• For the monostatic case both are equal and subscripts are not required. 


D. PHYSICAL OPTICS CURRENT 

The components of the surface current in terms of the local cartesian coordinate 
variables are from [Ref. 6] 


and 


Eg. cos Q” cos f Ej£ cos 
2 R s + Z a cos 6' 2 R s cos Q" + Z 0 


(3.37) 


_ ~ E ie . cos#' 1 'sin 0* Ej£ cos#"cos0* 
** ” I 2 R s Tz o cos#' 2 R s cos#' + Z 0 


(3.38) 


where R s is the surface resistivity. When R s = 0 (|T| = 1), the surface is a perfect 


electric conductor. As R s —» °° the surface becomes transparent (|r| —> 0). Defining the 


transverse magnetic (TM) and transverse electric (TE) reflection coefficients as 
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(3.39) 


r TM = r || = 
r TE = r i = 


~Z Q cos 9" 
2R S + Z 0 cos 6" 

ZlsL 

2 R s cos 6"+Z 0 


allows the surface current components to be expressed as 




E, r cosf^ , E if sinf r 

1 II T — 1 i 


/V/ 

COS# 


Jy = 


E ie - sin <j>" E l0 . cos f . 


■ r !- 


cos^' 


The spherical polar angles are known from the direction cosines 

sin ^ = V («') 2 + (^) 2 
cos^ = Vl-(«') 2 +(v") 2 


/ * ' 

( V 




tan (/)" = arctan 

l 

In the bistatic case, the incident angles are used in the above formulas. 


(3.40) 


(3.41) 

(3-42) 


(3.43) 


E. DIFFUSE CALCULATION 

Models may have a rough surface due to irregularities or material imperfection. 
Surface roughness introduces path differences relative to a smooth surface and thus 
contribute a diffuse scattering component to the RCS that is added to the coherent 
component. The random deviation of the surface results in diffuse scattering which 
yields a relatively angle dependent term to the RCS. First we assume that the surface 
deviation resemble a Gaussian shape as shown in Figure 7. The extent of the Gaussian 
shape is characterized by a correlation distance C (inmeters), which is the average 
distance at which deviations become uncorrelated. A large correlation distance implies a 
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slowly varying surface error, whereas a small correlation distance is associated with a 

2 

rapidly varying error. The variance of the irregularities is denoted as 8 (or standard 
deviation as 8) as illustrated in Figure 8. 



Figure 7. Gaussian Approximation of Rough Surface From [Ref. 1] 



— *■ * 


Figure 8. Correlation Distance, (a) Large Correlation Distance, 
(b) Small Correlation Distance From [Ref. 1] 


In the PO calculation, surface roughness is approximated by defining values for 
correlation distance and standard deviation. The average value of the RCS is 
approximately [Ref. 7] 


a 


ave 


4 ttA 1 
~ 


cos 2 6 e 


- 4 k'-s- 


4k 2 8 2 nC 2 f 
—expj- 



(3.44) 
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where A is the area of the facet and 6 is measured from the normal. This value 
simply added to the specular term 



F. THE SCATTERED FIELD OF A COLLECTION OF FACETS 

An arbitrary target can be approximated by a collection of triangular facets. 
Equation 3.16 must be evaluated for all facets, and the polarization factors and coordinate 
systems may differ for each facet. However, the integral over all facets will have the 
form of equation 3.17. To obtain the scattered field of the array of facets the following 
procedure is used: 

1. The incident field is specified in the global coordinates: angle of incidence 
and polarization (E i 0 ,E i( p). 

2. For each facet (indexed from c = l,2,...,N c ) the incident field is 

transformed to the local facet coordinate system. Also the angle of incidence is 
determined in the local facet coordinate system. 

3. The physical optics current is computed using the incident magnetic field, and 
the radiation integral evaluated over the facet surface to obtain the scattered field 
components (E s q*,E S 0 ') . 

4. The scattered field components are transformed back to the global coordinate 
system, thereby obtaining (E s Q,E s( p) for each facet. 

5. All of the N c transformed scattered field components are superimposed 

(vector summed) to obtain the total scattered field in the observation direction. 

6. The RCS is computed from the total scattered field from equation (2.3). 

A flow diagram of the procedure for the monostatic calculation is shown in Figure 9. For 
a bistatic calculation, the current does not have to be recomputed for each observation 
angle because the incident angle, and hence current, is fixed. [Ref. 8] 


25 



READ 

COORDINATES 
AND FACETS 


INITIALIZE PARAMETERS 



Figure 9. Flow Chart for Computing Monostatic RCS of an Array of Facets 



















IV. THE POFACETS GRAPHICAL USER INTERFACE 


POFACETS is a graphical user interface (GUI) radar cross section (RCS) 
prediction code based on the physical optics (PO) approximation and implemented in 
MATLAB. It is accurate when the target is large in terms of wavelength. Scattering 
objects are approximated by arrays of triangles (facets) and superposition is used to 
compute the total RCS of the model object as described in Chapter III. Multiple 
reflections are not included. The GUI consists of three modules: (1) Design Model , (2) 
Calculate Monostatic RCS, and (3) Calculate Bistatic RCS. 

Design Model provides a convenient way of modeling a complex object or loading 
and editing a previously designed model. The geometry of the model is defined using 
two sets of data. The first set defines the x,y,z coordinates of the vertices. The second 
set of data defines the component triangular facets of the object, their nodes, and 
illumination and resistivity characteristics. 

Calculate Monostatic RCS computes the total radar cross section of the target using 
a monostatic radar (co-located transmitter and receiver). The scattered field of each 
triangle is computed as if it was isolated and other triangles were not present. Multiple 
reflections, edge diffraction and surface waves are not included. Shadowing is 
approximately included by considering a facet to be completely illuminated or 
completely shadowed by the incident wave. Options for defining the roughness of the 
surface are available. 
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Calculate Bistatic RCS computes the total radar cross section of the target using a 
bistatic radar (transmitter and receiver at different locations). The incidence direction is 
fixed and the observation direction changes. 

To run the POFACETS GUI, the user must first run MATLAB and set the path or 
directory to where the POFACETS code is stored. This may be done in two ways. The 
first is by using the File Menu and selecting Set Path. This displays the Path Browser 
window and from here click on the Browse button. In the Browse for Folder window, 
select the folder where POFACETS is located then click the OK button. The second 
method is by using the change directory command at the MATLAB command prompt to 
where the code is stored such as cd. c:\pofacets. After setting the path or directory, typing 
pofacets at the command line will execute the GUI as shown in Figure 10. A brief 
explanation of the POFACETS GUI may be obtained if the user clicks on the Help 
button, which displays the POFACETS Help GUI as shown in Figure 11. 
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-J> RCS PREDICTION USING PHYSICAL OPTICS APPROXIMATION 




EkaoE. Garrido Jr. and David C. Jean 


Create a new target geometry or load; 
and edit a previously designed model. 


Design Model 


Compute fee monostatic radar crass 
section of fee modeLl 


Calculate Monostatic RCS 


Figure 10. POFACETS Main GUI 
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POFACETS Help 


IPOFACETS is a MATLAB graphical user interface (GOI) radar cross 


section prediction (JR.CS) code based on the physical optics (PO) 
approximation method. It is accurate for relatively high frequencies. 
Scattering objects are approximated by arrays of triangles (facets) 
and superposition is used to compute the total RCS of the model 
object 

Design Model provides a convenient way of modeling a complex object 
or to load and/or edit a previously designed and saved model The 
geometry of the model is defined using two sets of data. The first 
set defines the xyz coordinates of the vertices. The second set of 
data defines the component triangular meets of the object, its nodes, 
illumination and resistivity characteristics. 

Calculate Monostatic 1RC$ computes the total radar cross section of 
the target using a monostatic radar (co-located transmitter and receiver). 
The scattered field of each triangle is computed as if it were isolated 
and other triangles were not present. Multiple reflections, edge 
diffraction and surface waves, are not included Shadowing is approximate! 


Figure 11. POFACETS Help GUI 


A* DESIGNING AND PLOTTING A MODEL 

To illustrate the Design Model module of the GUI, we will design the AIM-9 
Sidewinder missile shown below in Figure 12 by modeling it with the set of coordinates 
in Table 1 and facet configuration in Table 2. 



Figure 12. AIM-9 Sidewinder from [Ref. 9] 
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Vertex # 

X (m) 

Y(m) 

Z(m) 



Y(m) 

Z(m) 

1 

0.000 

0.000 

0.00 

20 

-0.250 

0.065 

0.150 

2 

0.130 

0.000 

0.00 

21 

-0.250 

0.065 

0.000 


0.130 

0.130 

0.00 

22 

0.065 

0.000 

0.000 


0.000 

0.130 

0.00 

23 

0.065 

0.000 

0.200 


0.000 

0.000 

2.57 

24 

0.065 

-0.250 

0.150 

6 

0.130 

0.000 

2.57 

25 

0.065 

-0.250 

0.000 


0.130 

0.130 

2.57 

26" 

0.130 

0.065 

2.200 

8 

0.000 

0.130 

2.57 

27 

0.130 

0.065 

2.400 


0.065 

0.065 

2.87 

28 

0.380 

0.065 

2.200 

10 

0.130 

0.065 

0.00 

29 



2.200 

11 

0.130 

0.065 

020 

30 

0.130 

0.065 

2.400 

12 

0.380 

0.065 

015 

31 

0.065 

0.380 

2.200 

13 

0.380 

0.065 

0.00 

32 

0.000 

0.065 

2.200 

14 

0.065 






2.400 

15 

0.065 

0.130 

0.200 

34 

-0.380 

0.065 

2.200 

16 

0.065 

0.385 

0.150 

35 

0.065 

0.000 

2.200 

17 

0.065 

0.385 

0.000 

36 

0.065 

0.000 

2.400 

18 

0.000 

0.065 

0.000 

37 

0.065 

-0.380 

2.200 

19 

0.000 

0.065 

0.200 



Table 1. AIM-9 Sidewinder Coordinate Set 
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Facet # 

Node 1 

Node 2 

Node 3 

Illumination 

Normalized 

Resistivity 

1 

3 

2 

1 

1 

1 

2 

1 

4 

3 

i 

i 

3 

3 

7 

6 

i 

i 

4 

6 

2 

3 

i 

i 

5 

4 

8 

7 

1 

1 

6 

7 

3 

4 

i 

i 

7 

4 

1 

5 

i 

i 

8 

5 

8 

4 

i 

i 

9 

2 

6 

5 

1 

1 

10 

5 

1 

2 

1 

1 

11 

2 

9 

6~ 

1 

1 

12 

8 

9 

7 

i 

1 

13 

5 

9 

8 

1 

1 

14 

5 

6 

9 

1 

1 

15 

10 

11 

12 

0 

i 

16 

12 

13 

10 

0 

1 

17 

16 

15 

14 

0 

i 

18 

14 

17 

16 

0 

1 

19 


19 

18 

6 

i 

20 

18 

21 

20 

0 

1 

21 

22 

23 

24 

0 

1 

22 

24 

25 

22 

0 

1 

23 

26 

27 

28 

0 

i 

24 

29 

31 

30 

0 

i 

25 

32 

34 

33 

0 

i 

26 

35 

36 

37 

0 

1 


Table 2. AIM-9 Sidewinder Facet Configuration 
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To design a new object, the user clicks on the Design Model button of Figure 10. 
This displays the Design Model Coordinates GUI shown in Figure 13. 



Figure 13. Design Model Coordinates GUI 


The Number of Vertices is initially disabled but selecting File j New in the menu 
will enable this box to accept an entry. Since the Sidewinder example has 37 vertices, we 
enter this number in the appropriate box and the GUI generates the corresponding 
number of editable text boxes (all initialized to zero) to accommodate the Cartesian 
coordinates of the model as well as enable the Apply button to check all the individual 
coordinates later. Entering the coordinates correctly and clicking on the Apply button 
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displays the Coordinate Status GUI (Figure 14) and after clicking the OK button, enables 
the Design Facets button and the active GUI should now look like Figure 15. 


Coordinate Status 


mm 

Coordinates valid, press Design Facets button. 
OK | 


Figure 14. Coordinate Status GUI 


DESIGN MODEL Coordinates 


Fife .'E-EiE' 

Number of Vertices: J 37 1.j—7— 

' 4 VERTEX COORDINATES (m> SCaRK 1 1 

vF ..X . Y 2 v# X Y : W 

1.I. 1..0 M .o 1 , . 0} 

1.21 jv|._. -0.25 j ... 0.0651 C .0 

1 2 1 0.13 1 0 I 0 

Eli-E.' 0.035 C 0 f~ 0 

! 3 ,1 0.13 j 0.13 j 0 

I 23 lE. 0165 j.. 0/4..... .0.2 .S 

1 i/ I. 0:| .... 0.13 1 _0 * 

. .. 

1,1.1... 0 1 0 1 257 

1:25 m 0,065 1 *0.25 1 0 

■ El. 1., .0.13 j .0:4 257 

\.m 1_111; L_0.065 l ..22 

1 ? 1 0.13 J 013 J 257 

? fi i , 0 j 0.13' j 257 

rf* i... q- 131; i..." 0.06 = j:..::e4. 

Jlp 1.“.' 038 f~ , 0,065 r . 2.2 

1 S. |.....0,035 1_0065 j_287 ; 

. Oil 1 0.13 j 0065 1 0 

JW. 1 0.065; j~ ....0.13 t 2.2 

Tm 1 0055,. j ... ...0.13 . J 24 


JW f 0.055 1 0.38 1 22 

j -2 1 0.38 J 0.065 : • 1 0.15 

132 ].......0 J .0065 J__22 


1 33 1...... 0j4, ...1165 r .2.40. 

I G.0S5 j 0.13 1 0 

1 34 ) -038 ■ j 0065 1 22 

El 1...0.065 1_.0,13. L_1,2 

111: j.0,065 i . 0 i 22 

1 11; 1.. 0065. j.0.385. j.01.5 

1.23.4..0,065 J. 0 j. 24 

1 17! | 0,065 ., El . 0.385 1 0 

: : OE 1 o j 0.065 jE 0 

JIT j 0065 r -0.38 r 2.2 

1 19] j 0 1 0.065 1 0.2 

[20 J -025,. I 0065 | 0.15 

;.E Design Facets Apply : 

Close Help 


Figure 15. Design Model Coordinates GUI with Vertices of the Sidewinder 
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* DESIGN MODEL Facets 


















is enabled. The nodes and illumination values are then encoded using the data in Table 2. 
The illumination parameter specifies whether both sides of a facet are illuminated by the 
incident wave. If Hum = 0 both sides of the plate are illuminated (e.g., a fin or a wing); 
if Hum = 1 only the outer face is illuminated (e.g., a fuselage). The normalized resistivity 
R s !377 is used to define the resistivity characteristics of the model surface. For 
example, to conveniently define the entire missile skin to have an impedance of 377 
ohms, uniform resistivity is set to 1 and it will be noted that the resistivity values for all 
the facets are accordingly set in the Rs column. Pressing the Apply button displays the 
Facet Configuration Status GUI (Figure 17) if there are no errors and subsequently 
clicking the OK button enables the View Model button and the active GUI should now 
look like Figure 18. 


**# Facet Configuration Status 


Ena 



All nodes, illumination, end resistivity values valid. Click View Model 
Sutton. .■ ' •: 


OK 


Figure 17. Facet Configuration Status GUI 


To visualize the model comprised of its triangular components, the plot may be 
seen by clicking the View Model button, which displays the Model Geometry GUI as 
shown in Figure 19. 
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Figure 18. Design Model Facets GUI with Facet Configuration of the Sidewinder 



Figure 19. Model Geometry GUI 
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1. Help GUIs 

A number of Help GUIs are available to provide the user valuable assistance in 
working with the Design Model module. The three Help GUIs may be viewed by 
clicking on the Help buttons of the Design Model Coordinates, Design Model Facets, and 
the Model Geometry GUIs to display Figures 20, 21 and 22, respectively. 





Figure 20. Design Model Coordinates Help GUI 
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- # Design Model Facets Help 

j Design Model Facets conjures the geometry of the object by defining 


an array of triangular components or facets comprising the model The 
facets are described using a facet number, which by default is 
sequential. Each facet is defined by three nodes or vertices, and its 
illumination and resistivity characteristics. 

Nodes - Entries in the nodes (nodel, node2, and node3) represent the 
vertex numbers defining the facet The nodes comprising the facet 
should be ordered n a right-hand or counter clockwise sense where 
the normal to the facet is the direction of illumination. 

Bumination (Sum) - The ffiuminarion value is 0 if the facet is 
allowed to be illuminated from both sides and 1 if the meet will 
only be illuminated from one side. 


Resistivity (Rs) - The resistivity of the meet is normalized to the 
free space value of 37? ohms. If surface Rs = 0, then the facet is a 
perfect electric conductor (PEC); when Rs is infinity (inpractice, 
several thousand) the facet is transparent 




Close 


Figure 21. Design Model Facets Help GUI 


* Model Geometry Help 


The plot of the Model Geometry is a scaled version of the design 


defined by the coordinates and the facets configuration. 

Label Vertices - is a function that plots the model geometry showing 
the vertex numbers as defined in the Design Target Model-Coordinates 
GDI 

Labe! Faces - is a function that plots the model geometry showing the 
facet or triangle number as defined by the nodes in the Design Target 
Model - Facets. 

Grid On/Off- is a function that plots or clears the axis area of the 
model geometry figure of grid fines to enhance visualization of the 
model. 

Print Files - this function prints a hard copy of the c o ordinate s/£ac et s 
data and/or the model geometry as selected 

Close - doses the GDI figure. 




l'}A' 


vVy; 




wv. 


.Close 


Figure 22. Model Geometry Help GUI 
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2 . 


Status GUIs 


Displaying several Status GUIs have been incorporated in the code to draw user 
attention to specific part of the GUI requiring action. For the Design Model Coordinates 
GUI, an error in the Number of Vertices entered will display the status shown in Figure 
--*■ Likewise, since the Design Model Coordinates GUI can initially accommodate a 
maximum of 40 vertices only, exceeding this limit will display Figure 24. Meanwhile, 
encoding an invalid Scale value will show' Figure 25. 


I Number of Vertices Status 

sssi 

Please check the numb 

. OK j 

er of vertices entered! u 


Figure 23. Number of Vertices Status GUI 


Maximum Number of Verices 


as m 


A Design can handle up to 40 vertices only! 


: OK 


Figure 24. Maximum Number of Vertices Status GUI 


I Scale Status 

mmmw 

O 

: ‘lease check scale of vertices entered. i 

; OK | 


Figure 25. Scale Status GUI 
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The user is advised of erroneous entries on the coordinates with the display of 
Figures 26a to c, while a valid set of coordinates is ascertained by Figure 26d. 




(d) 


Figure 26. Coordinate Status GUIs, (a) Invalid Coordinate, (b) Uncorrected Erroneous 
Coordinate, (c) Duplicate Set of Coordinates, (d) Valid Set of Coordinates 











For the Design Model Facets GUI, an error in the Number of Facets entered will 
display Figure 27. Likewise, since the GUI can initially accommodate a maximum of 40 
facets only, exceeding this limit will display Figure 28. Meanwhile, encoding an invalid 
resistivity value for a particular facet will show Figure 29. 



Figure 29. Uniform Resistivity Status GUI 


The user is advised of erroneous entries in the configuration of the facet(s) with 
messages shown in Figures 30a to f and encoding valid parameters is ascertained by 
Figure 30g. 
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Figure 30. 



(a) 


Node Status 


0 


I ■■■ 3Ifl 

Duplicate nodes for the same facet number! 
OK i 


(b) 


Node Status 


□m 



Error in node(s), illumination or resistivity entered! Correct entries 
where 9993 is displayed. 


OK 




(C) 



(d) 


Facet Configuration Status GUIs, (a) Invalid Node, (b) Duplicate Nodes, 
(c) Unconnected Erroneous Node, (d) Duplicate Set of Nodes 
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(§) 

Figure 30 (continued). Facet Configuration Status GUIs, (e) Invalid Illumination Value, 
(f) Invalid Facet Resistivity, (g) Valid Facet Configuration 


3. Model Geometry GUI Functionalities 

The Model Geometry GUI shown in Figure 19 provides some functionalities to 
enhance user visualization of the model object and to verify correctness of node and facet 
designation. These include plotting the vertex/facet numbers with the object or the grid 
lines by selecting any or all of the following checkboxes: Label Vertices, Label Faces, 
and Grid, respectively. Selecting the Coordinates/Facets and Model Geometry 
checkboxes and clicking on the Print button will generate a hard copy of the coordinates 
and facets data. 
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4. Saving the Model Files 

The model’s file is composed of two sets of data, coordinate set and facet 
configuration. However, to save a model only the filename is required. To save the 
Sidewinder example, click on File | Save As... from the menu and this displays the New 
File GUI shown in Figure 31. We enter the filename Sidewinder and click on the Save 
button. The program saves the files and confirms the action to the user by displaying the 
File Saved GUI shown in Figure 32. Inadvertently specifying a filename that already 
exists in the list of models will prompt the user to overwrite said file with the Replace 
File GUI shown in Figure 33. 






£ 

.. 

is?,"": 


BOX 

1 

WAIL 

—I 

PLATE 


BRICK 


CHUCK_ . _ 

zl 

File name: 

Save | 

Sidewinder] 

Cancel 


Figure 31. New File GUT 



Figure 32. File Saved GUI 
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c:^o1aeels^odels^LONG_BOXaIready exists. 
Do you want to replace It? \ 



Figure 33. Replace File GUI 


5. Editing an Existing Model 

To edit an existing model, its set of data needs to be retrieved and displayed 
before editing or scaling can proceed. This is accomplished by clicking on the File | 
Open menu of the Design Model Coordinates GUI (Figure 10) which displays the Open 
Model File GUI shown in Figure 34. Suppose we would like to edit the VTAIL model. 
We select it among the listed file names and click the Open button. This action performs 
a series of operations that will display its data set and the model plot as shown in Figure 
35. 



Figure 34. Open Model File GUI 
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Figure 35. Displayed VIAIL File 


At this point, the user may click on any of the three overlapping GUIs to make it 
active and edit the pertinent data. The different Status GUIs mentioned earlier will guide 
the user throughout the use of the interface. 

B. CALCULATING THE MONOSTATIC RCS 
1. Smooth PEC Model 

To compute the Monostatic RCS of a model object, the user clicks on the 
Calculate Monostatic RCS button of the POFACETS Main GUI (Figure 10). This clears 
the screen of other GUIs (except the Main GUI) and brings up the Calculate Monostatic 
RCS GUI with default values and parameters shown in Figure 36. 
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Figure 36. Calculate Monostatic RCS GUI 


For this case, we will calculate the monostatic RCS of a closed box model with a 
smooth surface, measuring 100 meters on each side and which is made of perfect electric 
conductor (PEC) material. Since this is a closed form object, only the outer face will be 
illuminated and all the facets will have a uniform illumination of 1. We further assume 
that the incident wave is 0-polarized (TM- z ), the frequency of interest 15 MHz and we 
would like the RCS measured at o — 0 plane while varying 6 from 0 to 360 degrees to 
get a view of the model RCS around its sides. 

Setting up the GUI with the parameters mentioned should result in the Calculate 
Monostatic RCS GUI that looks like Figure 37 and the RCS plot shown in Figure 38. 
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Two computational parameters that may be varied are the Taylor Series Region and the 
Number of Terms in the Taylor Series. See [Ref. 5] for details on how these parameters 
affect the accuracy and computation time. In general, the default values give good 
results with relatively fast run times. Also, a two-dimensional contour can be generated 
by setting 6 values from 0 to 90, the <b values from 0 to 360, increment angle to 2 
degrees, and the resulting computation of the monostatic RCS is shown in Figure 39. 


¥1 CALCULATE MONOSTATiC RCS 


Model Name: 



#i 


BRICK 

CHUCK 

DART 

; : : : : 

LONG BOX 


PLATE 


Sidewinder 


UTA2 

w 

VTAIL 

.. .3 


Lead File 


Surface Roughness 

Coirdaiion Dist.(m).j 0 
Standard Dew. (m): I 0 


Theta 


Phi 


Starting Angie >■ j 0 

deg 

1 0 

deg 

Ending .Angle | 360 

"deg T 

1 0 

deg 

increment Angle 2 

; de g' 

i 

deg 


5 Computational Parameters 

Taylor Series 

Length of Region: J C.Cd 
N umbercfTerms: "S 5 


Incident Polarization 
j Theta (TM-z) zl 


Frequency j 0.015 

.■.*. 4 ;:;.:. . GHz ' 


Calculate RCS 


Print 


Close 



Figure 37. Calculate Monostatic RCS GUI for Smooth PEC Box Model 


49 

















Figure 38. Monostatic RCS of Smooth PEC Box Model 



* MONOSTATIC RCS 


Figure 39. Contour Plot of Smooth PEC Box Model 
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RCS-theta is specified as Gee (# polarized incident wave and # polarized receiver). 
This is the co-polarized component. RCS-phi is o^e (0 polarized incident wave and <p 
polarized receiver). The latter is the cross-polarized component. In this particular 
example, it happens to be at a level below the specific threshold in the code. 

2. Model with Surface Impedance 

An object coated with a radar absorbing material or having a skin made of 
resistive material may be approximated by defining its surface to have an impedance 
value. Equations 3.37 and 3.38 describes the components of the surface current at the 
target which is inversely proportional to the resistivity of the surface. Using the box 
model as an example and defining its surface to have an impedance three times the value 
of free space (i.e. uniform resistivity is set to 3), we see that the monostatic RCS 
calculation (Figure 40) is 15dB down in the specular return and from 15 to 20dB at the 
other various angles compared to the PEC model shown in Figure 38. 

3. Rough Surface Model 

Meanwhile, to observe the effect of the diffused scattering component of a rough 
surface, the box model is defined to have a correlation distance of 100 meters (or 
equivalently 5A.) and a standard deviation about the mean value of 2 meters (that is, 
0. 1A for a frequency of 15 MHz). Without varying the other parameters, the monostatic 
RCS plot in Figure 41 shows a decline of about 8 dB in the specular direction and an 
increase in the RCS at all the other angles due to the diffuse scattering effect compared to 
the PEC model results in Figure 38. 
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Figure 40. Monostatic RCS of Smooth Box Model with Surface Impedance 



Figure 41. Monostatic RCS of Box Model with Rough Surface 
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Calculate Monostatic RCS GUIs 


Several GUIs were designed to enable the user to pinpoint specific parameters 
requiring attention. Errors in encoding values for the theta and phi angles display Figures 
42a to c, while setting the same starting and ending angles advise the user that a theta or 
phi-cut is set as shown in Figure 43. 



(c) 

Figure 42. Angle Status GUIs, (a) Invalid Theta Starting Angle, (b) Invalid Theta Ending 

Angle, (c) Invalid Theta Increment Angle 
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Figure 43. Theta-Cui Set GUI 


Other GUIs include notifying the user of the need to correct entries for the 
Frequency, Wavelength, and the Length of the Taylor Series Region shown in Figures 44, 
45, and 46, respectively. 



Figure 45. Wavelength Status GUI 



Figure 46. Taylor Series Region GUI 
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In addition, some GUIs have been pre-defined to warn the user if the parameter 
entered will take some time to generate the results as in Figure 47 or suggest a range of 
values such as Figures 48 and 49. 


Taylor Series Terms 


:.lslx| 



Monosiatic RCS computation for number of terms greater then 10 may take 
sometsime. : '■? : A- V' 


OK 


Figure 47. Taylor Series Terms GUI 


U] Correlation Distance Status 





T 

* 


Suggest that this value be smaller than the wavelength. 

OK 


Figure 48. Correlation Distance Status GUI 


Standard Deviation Status 


nsm 



A\ Value entered Is valid bii extremely high. 


K i: OK 


Figure 49. Standard Deviation Status GUI 


It will be evident to the user that the Calculate RCS button is initially disabled. 

Selecting a model from the list box and clicking the Load File button will confirm to the 

user the model selected and give a reminder to set all the parameters as shown in Figure 

50. Clicking the OK button will enable the Calculate RCS button. Finally, clicking the 

Calculate RCS button informs the user that the GUI is calculating the monostatic RCS of 
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the model selected as shown in Figure 51 below. The Calculate Monostatic RCS Help 
GUI showm in Figure 52 can also provide some valuable information. 




Figure 51. Computing Monostatic RCS GUT 


Calculate Monostatic RCS Help 


f Calculate Monostatic RCS computes the total radar cross section of a 


■ model using a monostatic radar (collocated transmitter and receiver). 
Scattering objects are approximated by arrays of triangles (facets) and 
trie scattered field of each triangle is confuted as if it were isolated 
and other triangles were not present. Multiple reflections, edge 
difraction and surface waves are not included. Shadowing is 
approximately included by considering a facet to be completely 
: illuminated or completely shadowed by trie incident wave. 


Start and Ending Angles - The illuminated surface of a model is described 
by defining the starting and ending angles (theta and phi values) in degrees 
: of the model with respect to trie monostatic radar. Angles are measured in 
degrees and entering a phi start and phi stop with equal values results in a 
phi-cut. Similarly, entering a theta start and theta stop with equal values 
: results in a theta-cut. 

: Increment - defines the number of degrees (or fraction of a degree) increase 
in step for the start and ending angles to be used in calculating the mono static 
radar cross section of the model.__ 


Figure 52. Calculate Monostatic RCS Help GUI 
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C. CALCULATING THE BISTATIC RCS 

The Calculate Bistatic RCS GUI (Figure 53) is very similar to the Calculate 
Monostatic RCS GUI (Figure 36) except for the 6 and 6 incident angles that define the 
transmitter location with respect to the object. 



Figure 53. Calculate Bistatic RCS GUI 


For an incidence angle of 45 degrees in the $ = 0° plane and the observation 
direction varied from 0 to 360 degrees in the <f> = 0° plane, using a frequency of 0.2 GHz 
and the number of Taylor Series terms set at 3, the bistatic RCS of the box model is 
shown in Figure 54. The lobe at 135 degrees is the specular reflection from the top. The 

forward scatter is also evident. Note that the forward scatter is higher than the highest 
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monostatic backscatter lobe. The accompanying ray diagram in Figure 53 illustrates the 
directions of the expected large RCS lobes. 















Status GUIs similar to that of the Calculate Monostatic RCS module were 
designed to enable the user to easily correct specific parameters. The Calculate Bistatic 
RCS, however has an additional GUI for the incidence angle as shown in Figure 56. 
Other helpful tips may be found in Calculate Bistatic RCS Help GUI shown in Figure 57. 


# Incident Angle Status 


m 


© 


Enter an incident angle between -360 and +360 degrees. 


OK 


Figure 56. Incident Angle Status GUI 


BISTATIC RCS Help 


Calculate Bistatic RCS computes the total radar cross section of the target P 


using a bistatic radar (transmitter and receiver at different locations). 

Incidence Angle - The incidence direction is fixed and is defined by a theta 
and phi angle in degrees. 

Observation Angles - "WHe the incidence angle is fixed, the observation 
direction changes. Theta/Phi Start and Stop Angles are measured in degrees 
and entering a phi start and phi stop with equal values results in a phi-cut. 
Similarly, entering a theta start and theta stop with equal values results in a 
theta-cut. 

Increment - defines the number of degrees (or fraction of a degree) increase 
in step for the thetafchi start and stop angles to be used in calculating the 
bistatic radar cross section of the model. 

Surface Roughness - A model may be designed to have a rough surface by 
inputting values for correlation distance and standard deviatioa These 
parameters contribute a disuse scattering component of the RCS and i s adde 


Close 


Figure 57. Calculate Bistatic RCS Help GUI 
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V. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

The main purpose of this work was to develop a Graphical User Interface (GUI) 
that will provide a convenient way of modeling complex objects using triangular facets, 
provide options to define the surface roughness of the scattering object, and calculate the 
monostatic or bistatic radar cross section over a range of frequencies or wavelengths of 
interest. Care was taken to design the GUI to make it appealing and easy to use and 
several validation techniques were incorporated to make it robust and prevent the user 
from encoding erroneous entries. Help windows and messages guide the user in making 
efficient use of the GUI. The GUI was implemented in MATLAB utilizing its scientific 
computing features and several built-in plotting functions to visualize the model 
geometry and results of the RCS calculation. Implementing the code as a user interface 
is the best way to avoid errors that may have significant consequences related to the use 
of this tool for the desired purpose. 

The GUI can be used by both students and electromagnetic professionals to 
quickly determine the RCS performance of their model by varying the computational 
parameters to suit their design requirements. 

B. RECOMMENDATIONS 

The GUI developed was based on the basic physical optics approximation. An 

important area for future work is to be able to upgrade this GUI into one that incorporates 

other scattering mechanisms, primarily the effects of multiple reflections. Multiple 

reflections complicate the PO approximation because the incident field on a patch is not 

only determined by the direct field from the source, but also the field scattered by other 
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facets in the direction of the triangle under consideration. Numerical integrations will be 
required that substantially increase the computation time. 
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APPENDIX A. MATLAB M-FILES 


The following are the description of the Matlab script files contained in the 
POFACETS Graphical User Interface: 

1. BPlot.m 

This file generates either a theta of phi-cut bistatic RCS calculation plot of a 
selected model. 

2. BPlotl.m 

This file generates a contour plot of the bistatic RCS calculation of a selected 

model. 

3. BMessage.m 

This file displays the name of the model selected in the Calculate Bistatic RCS 
GUI and enables the Calculate RCS button. 

4. G.m 

This file calculates the recursive Green’s function using the two parameters (n and 
y) passed. 

5. MPIot.m 

This file generates either a theta of phi-cut monostatic RCS calculation plot of a 
selected model. 

6. MPlotl.m 

This file generates a contour plot of the monostatic RCS calculation of a selected 

model. 
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7. 


MaxPoints.m 


This file displays a message that the GUI design can accommodate only a 
maximum of 40 vertices only. 

8. Repl_File.m 

This file displays a GUI prompting the user to replace the data of an existing 

model. 

9. Save_As.m 

This file displays a GUI with a list of the existing models and prompts the user to 
enter a filename. 

10. SvAs.m 

This file implements the functionality of the Save_As GUI by displaying the 
filename selected from the list and/or checking if the filename entered already exists. 

11. VertS tatus.m 

This file displays a GUI that prompts the user to check the number of vertices if 
the value entered is erroneous. 

12. bistatic.m 

This file displays the Calculate Bistatic RCS GUI. 

13. bistat.m 

This file implements the functionalities of the Calculate Bistatic RCS GUI by 
checking all the parameters entered, displaying GUIs with appropriate messages, 
calculating the bistatic RCS of the selected model, and generating the plot. 
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14. design.m 

This file displays the Design Model Coordinates GUI for encoding the 
coordinates of the model being designed or edited. 

15. dsgn.tn 

This file implements the functionalities of the Design Model Coordinates GUI by 
checking the individual coordinates entered and the whole coordinate set. 

16. dsgnFN.m 

This file clears previous design GUIs and displays a new Design Model 
Coordinates GUI. 

17. dsgnFO.m 

This file displays the Open Model File GUI showing the list of existing models 
which the user may select to display the data set and plot of the geometry. 

18. facets.m 

This file displays the Design Model Facets GUI for defining the configuration of 
the component triangular facets of the model. 

19. fact.m 

This file calculates the factorial of a number n given the parameter. 

20. fcets.m 

This file implements the functionalities of the Design Model Facets GUI by 
checking the number of facets, individual nodes and set of nodes, illumination and 
resistivity values. 
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21. hlpbistattxt.m 

This file contains the text that is displayed in the Calculate Bistatic RCS Help 

GUI. 

22. hlpcoordtxt.m 

This file contains the text that is displayed in the Design Model Coordinates Help 

GUI. 

23. hlpfacettxt.m 

This file contains the text that is displayed in the Design Model Facets Help GUI. 

24. hlpmaintxt.m 

This file contains the text that is displayed in the POFACETS Main GUI. 

25. hlpmodeltxt.m 

This file contains the text that is displayed in the Design Model Geometry Help 

GUI. 

26. hlpmonotxt.m 

This file contains the text that is displayed in the Calculate Monostatic RCS Help 

GUI. 

27. hlpbistabm 

This file displays the Calculate Bistatic RCS Help GUI. 

28. hlpdsgncoord.m 

This file displays the Design Model Coordinates Help GUI. 

29. hlpdsgnfacet.m 

This file displays the Design Model Facets Help GUI. 
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30. hlpmain.m 

This file displays the POFACETS Main Help GUI. 

31. hlpmodel.m 

This file contains the text that is displayed in the Design Model Facets Help GUI. 

32. hlpmono.m 

This file displays the Calculate Monostatic RCS Help GUI. 

33. message.m 

This file displays the name of the model selected in the Calculate Monostatic RCS 
GUI and enables the Calculate RCS button. 

34. mod.m 

This file implements the functionalities of the Model Geometry GUI by plotting 
the model, labeling its vertices and facets, displaying the grid, and printing the data set 
and model plot. 

35. model.m 

This file displays the Model Geometry GUI. 

36. mono.m 

This file implements the functionalities of the Calculate Monostatic RCS GUI by 
checking all the parameters entered, displaying GUIs with appropriate messages, 
calculating the monostatic RCS of the selected model, and generating the plot. 

37. monostatic.m 

This file displays the Calculate Monostatic RCS GUI. 
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38. mplt.m 

This file prints the Monostatic or Bistatic RCS plot. 

39. msgComputing.m 

This file displays a message that the GUI is calculating either the monostatic or 
bistatic RCS of the model selected. 

40. pofacets.m 

This file displays the POFACETS Main GUI where the user may choose to run 
any of the three modules - Design Model, Calculate Monostatic RCS and Calculate 
Bistatic RCS. 
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